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£\j , We construct a theory for the 1+1-dimensional Brownian motion in a viscous medium, which is 

(i) consistent with Einstein's theory of special relativity, and (ii) reduces to the standard Brownian 
motion in the Newtonian limit case. In the first part of this work the classical Langevin equations 
of motion, governing the nonrelativistic dynamics of a free Brownian particle in the presence of a 
heat bath (white noise), are generalized in the framework of special relativity. Subsequently, the 
0^ ' corresponding relativistic Langevin equations are discussed in the context of the generalized Ito 

(pre-point discretization rule) vs. the Stratonovich (mid-point discretization rule) dilemma: It is 
found that the relativistic Langevin equation in the Hanggi-Klimontovich interpretation (with the 
post-point discretization rule) is the only one that yields agreement with the relativistic Maxwell 
distribution. Numerical results for the relativistic Langevin equation of a free Brownian particle are 
■ presented. 
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For almost one-hundred years, Einstein's theory of special relativity 0,0 is serving as the foundation of our most 
successful physical standard models (apart from gravity). The most prominent and, probably, also the most important 
feature of this theory is the absolute character of the speed of light c, representing an unsurmountable barrier for 
' the velocity of any (macroscopic) physical process. Due to the great experimental success of the original theory, 
almost all other physical theories have successfully been adapted to the framework of special relativity over the past 
decades. Surprisingly, however, the scientific literature provides relatively few publications on the subj ect of relativistic 
Brownian motions (classical references are 0,0,0, more recent contributions include [il OHMS EI El EI El)- 
H ' Brownian particles are physical objects (e.g., dust grains), which move randomly through a surrounding medium 
j— I ' (heat bath). Their stochastic motions are caused by permanent collisions with much lighter constituents of the 
heat bath (e.g., molecules of a liquid). The classical theory of Brownian motion or nonrelativistic diffusion theory, 
respectively, was developed by Albert Einstein an d Marian von Smoluchowski El • Since the beginning of the 
last century, when their seminal papers were published, the classical theory has been investigated and generalized by 
a large number of physicists [p| El El EI E3 an d mathematicians [2p |2^, E|. The intense research led, among 
others, to different mathematical representations of the Brownian motion dynamics (Langevin equations, Fokker- 
Planck equations, etc. [l8l liH . l20 | ) . to the notion of Wiener processes j2^|, and to new techniques for solving partial 
differential equations (Feynman-Kac formula, etc. [22u23| ). 

With regard to special relativity, standard Brownian motion faces the problem that it permits velocity jumps Av, 
that exceed the speed of light c (see also Schay 3J). This is due to the fact that in the nonrelativistic theory the 
velocity increments Av have a Gaussian distribution, which always assigns a nonvanishing (though small) probability 
to events Av > c. This problem is also reflected by the Maxwell distribution, which represents the stationary velocity 
distribution for an ensemble of free Brownian particles and permits absolute velocity values v > c |20| . 

The first relativistically consistent generalization of Maxwell's velocity distribution was introduced by Jiittner [24| 
in 1911. Starting from an extremum principle for the entropy, he obtained the probability distribution function of 
the relativistic ideal Boltzmann gas (see Eq. I|67[) below). In principle, however, J iittner's approach made no contact 
with the theory of Brownian motion. Fifty years after Jiittner's work, Schay Q performed the first comprehensive 
mathematical investigation of relativistic diffusion processes, based on Lorentz-invariant transition probabilities. On 
the mathematical side, Schay's analysis was complemented by Hakim || and Dudley 0|, who studied in detail the 
properties of Lorentz-invariant Markov processes in relativistic phase space. After forty more years, Franchi and Le 
Jan |13| have recently presented an extension of Dudley's work to general relativity. In particular, these authors 
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discuss relativistic diffusions in the presence of a Schwarzschild metric 25]. Hence, over the past one-hundred years 
there has been steady (though relatively slow) progress in the mathematical analysis of relativistic diffusion processes. 

By contrast, one finds in the physical literature only very few publications which directly address the topic of the 
relativistic Brownian motion (despite the fact that relativistic kinetic theory is fairly well established for more than 
thirty years |26l 1271 1281 29] ) . Among the few exceptions are the papers by Boyer |g, |9j] and Ben-Ya'acov |(| , who have 
studied the interaction between two energy-level particles and electromagnetic radiation in thermal equilibrium, the 
latter acting as a heat bath. In contrast to their specific microscopic model, we shall adopt a more coarse-grained 
point of view here by assuming that the heat bath is sufficiently well described by macroscopic friction and diffusion 
coefficients. 

Generally, the objective of the present paper can be summarized as follows: We would like to discuss how one can 
construct, in a physically straightforward manner, a relativistic theory of Brownian motion for particles moving in a 
homogeneous, viscous medium. For this purpose it is sufficient to concentrate on the case of 1+1-dimcnsions (gener- 
alizations to the 1+3-dimensions are straightforward and will be discussed separately in a forthcoming contribution). 
As a starting point we choose the nonrelativistic Langevin equations of the free Brownian particle. In Sec. ITTI these 
equations will be generalized such that they comply with special relativity. As we shall see in Sec. Mil due to multi- 
plicative noise for the momentum degree of freedom, the resulting relativistic Langevin equations are not sufficient in 
order to uniquely determine the corresponding Fokker-Planck equation (generalized Ito-Stratonovich dilemma). Fur- 
thermore, it it is shown that the stationary solution of a particular form for the relativistic Fokker-Planck equation 
coincides with Juttner's relativistic Maxwell distribution (Sec. MI B 31) . Finally, we also discuss numerical results for 
the mean square displacement in Sec. IIVI 

It might be worthwhile to emphasize that the systematic Langevin approach pursued below is methodically different 
from those in Refs. 0, 0, H B II II H ES IH 111 E3 and also from the kinetic theory approach [H H3, |H H| . 
It is therefore satisfactory that our findings are apparently consistent with rigorous mathematical results, obtained 
by Schay and Dudley Q for the case of free relativistic diffusion. Moreover, it will become clear in Sec. IIVI that 
numerical simulations of the relativistic Langevin equations constitute a very useful tool for the numerical investigation 
of relativistic diffusion processes, provided that the discretization rule is carefully chosen. 

II. LANGEVIN DYNAMICS 

First the main properties of the nonrelativistic Langevin equations for free Brownian particles are briefly summarized 
(Sec. Ill A|) . Subsequently, we construct generalized Lorentz-covariant Langevin equations (Sec. Ill Bp . Finally, the 
covariant Langevin equations will be rewritten in laboratory coordinates (Sec. Ill Cj l. 

The following notations will be used throughout the paper: Since we confine ourselves to the 1+1-dimensional 
case, upper and lower Greek indices a,j3, ... can take values 0,1, where '0' refers to the time component. The 
1+1-dimensional Minkowski metric tensor with respect to Cartesian coordinates is taken as 

(^/ 3 ) = (r? a/3 )=diag(-l,l). 

Moreover, Einstein's summation convention is invoked throughout. 

A. Physical foundations 

Consider the nonrelativistic one-dimensional motion of a Brownian particle with mass m that is surrounded by 
a heat bath (e.g., small liquid particles). In the Langevin approach the nonrelativistic dynamics of the Brownian 
particle is described by the stochastic dynamical equations (see, e.g., [2(| Chap. IX) 

v(t) (la) 

-umv(t) + L(t), (lb) 

where v is the viscous friction coefficient. The Langevin force L(t) is characterized by 

(£(«)> = 0, (L(t)L(s))=2D6(t-s), (2) 

with all higher cumulants being zero (Gaussian white noise), and D being constant. More general models may include 
velocity-dependent parameters v and D (see, e.g., ^| H3, yjl H2 ) , but we shall restrict ourselves to the simplest case 
here. It is worthwhile to summarize the physical assumptions, implicitly underlying Eqs. Q: 



da; . 
dt 
dv 
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• the heat bath is homogeneous, 

• stochastic impacts between the Brownian particle and the constituents of the heat bath occur virtually uncor- 
rected, 

• on the macroscopic level, the interaction between Brownian particle and heat bath is sufficiently well described 
by the constant viscous friction coefficient v and the white noise force L, 

• Eqs. QJ hold in the rest frame E of the heat bath (corresponding to the specific inertial system, in which the 
average velocity of the heat bath vanishes for all times t) . 

In the following So will also be referred to as laboratory frame. 
In the mathematical literature, Eq. Ijlbf) is usually written as 

d [mv(t)} = -vmv{t)dt + dW(t), (3a) 

where W(t) is a one-dimensional Wiener process fl9l I22I |23| . i.e., the density of the increments 

w{t) = dW(t) = W(t + dt) - W{t) (3b) 

is given by 

w{tf 



P>(i)] 



V4ttL> dt 



exp 



ADdt 



(3c) 



Here the abbreviation w = dW has been introduced to simplify the notation in subsequent formulae. From Eq. i|3c|) 
one finds in agreement with @ 

Wi ,)=o, m<>»m> ={° D<U> (4) 

Depending on which notation is more convenient for the current purpose, we shall use below cither the physical 
formulation (JTJ or the mathematical formulation J3J). The two formulations can be connected by (formally) setting 

w(t) = dW(t) = L{t)dt. (5) 



B. Relativistic generalization 



It is well-known that in inertial coordinate systems, which are comoving with a particle at a given moment t, the 
relativistic equations must reduce to the nonrelativistic Newtonian equations (see, e.g., [23 Chap. 2.3). Therefore 
our strategy is as follows: Starting from the Langevin equations \Q or l|3al) . respectively, we construct in the first step 
the nonrelativistic equations of motion with respect to a coordinate frame £*, comoving with the Brownian particle 
at a given moment t. In the second step, the general form of the covariant relativistic equations motions are found 
by applying a Lorentz transformation to the nonrelativistic equations that have been obtained for S*. 

It is useful to begin by considering the deterministic (noise-free) limit case, corresponding to a pure damping of the 
particle's motion. This will be done Sec. Ill B II Subsequently, the stochastic force is separately treated in Sec. Ill B 21 



1. Viscous friction 

Setting the stochastic force term to zero (corresponding to a vanishing temperature of the heat bath), the nonrela- 
tivistic Eq. I|lbjl simplifies to 

dv . . . . . . 

m—(t) — —umvvb). (6) 
dt 

The energy of the Brownian particle is purely kinetic, 

E(t) = TO®-, (7) 
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and, by virtue of QB}, its time derivative is given by 

dE dv 2 

— — = mv— = -emu . (8) 
dt dt w 

As stated above, in the nonrelativistic theory the last three equations are assumed to hold in the rest frame Eo of the 
heat bath. Now consider another inertial coordinate system E„, in which the Brownian particle is temporarily at rest 
at time t or t* = t*(t), respectively, where i* denotes the £*-time coordinate. That is, in E* we have at time t 

v*(t) = v*(t*(t)) =0. (9) 

(Conventionally, we use throughout the lax notation g*(t) = g*(t*(t)), where g* is originally a function of i*.) With 
respect to the comoving frame £*, the heat bath will, in general, have a non- vanishing (average) velocity V*. Then, 
using a Galilean transformation we find that Eq. 10 in ^-coordinates at time t reads as follows 

m T^(*) = -vm[v*(t)-V*] i vmV*. (10a) 

OX* 

Similarly, in S^-coordinates Eq. JSJ is given by 

^ (t) =-vmv* (t) [v m (t) - V*] = 0. (10b) 
dt* 

Note that in the nonrelativistic (Newtonian) theory the left equalities in Eqs. I|10|) are valid for arbitrary time t. 
By contrast, in the relativistic theory these equations are exact at time t only if £» is comoving at time t. In the 
latter case, we can use Eqs. (|10|) to construct relativistically covariant equations of motion. Introducing, as usual, 
the proper time r by the definition 



Hr :d/,/l-^=d M /i..--^ ,11) 



and combining momentum = mv* and energy into a 1+1-vector (p") = (p°,p*) = (S*/c,p*), we can rewrite Eqs. 
(|10fl in the covariant form 

^ - /**) (/") = -my (0, «. - K) • (12) 

Let (u") and (U?) denote the 1+1- velocity components of Brownian particle and heat bath, respectively. Now it is 
important to realize that the covariant force vector f a can not be simply proportional to the 1+1-velocity difference, 

/: ^ -mi/« - K), (13) 

since, in general, at time t in £„ 

u1-U°= , C - - C ic- - C ^0. (14) 

y/l - vl/c 2 v 71 ~ K 2 /c 2 v/1 - K 2 /c 2 

However, we can write f™ in a manifestly covariant form, if we introduce the friction tensor 



which allows us to rewrite l|12l) as 



dn a 

-mvSptf-US). (16) 



dr 



This equation is manifestly Lorentz-invariant and we drop the asterisk from now on, while keeping in mind that the 
diagonal form of the friction tensor l|15f) is linked to the rest frame £* of the Brownian particle. In this respect the 
friction tensor is very similar to the pressure tensor, as known from the relativistic hydrodynamics of perfect fluids 
(see, e.g., [2|| Chap. 2.10). This analogy yields immediately the following representation 



(17) 
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It is now interesting to consider Eq. (|16fl in the laboratory frame T,q, denned above as the rest frame of the heat 
bath. There we have 

(^) = (c,0), (u") = (7c,7»), dr=-, 7^ 71 1 2 / 2 (18) 

Combining Ijl6(l . (|17|l . and l|18J) we find that the relativistic equations of motion in Eo are given by 

dp mv 

-r = -v , = (19a) 

di ^1 -v 2 /c 2 

dE = v mv 2 

dt ^l~v 2 /c 2 

Upon comparing (|19a|) with JBJ and l|19b|) with JSJ, one readily observes that the relativistic equations l|19|l do indeed 
reduce to the known Newtonian laws in the limit case v 2 /c 2 <C 1. 
Using the relativistic definitions 

E = -fmc 2 , p — jmv, (20) 

Eqs. ifHfll can also be rewritten as 

| = (21a) 

d£ u 2 . 

■ = -vpv = -vE— (21b) 
di cr 

In fact, only one of the two Eqs. (|19fl or (|21|) . respectively, must be solved, due to the fixed relation between relativistic 
energy and momentum: 

M « = -£yi P 2 = -™V =* g(t)= , mc2 2 / 2 . (22) 

V 1 — w^/c^ 

The solution of (|21a|l reads 

p(i) =p exp(-W), p(0)=p o , (23) 

and, by using (|2U|I . one thus obtains for the velocity of the particle in the laboratory frame So (rest frame of the heat 
bath) 



v(t) = v 



1-1)^ + 1 



-1/2 



(24) 



Figure depicts a semi-logarithmic representation of the velocity v(t) for different values of the initial velocity Vq. 
As one can see in the diagram, at high velocities \v\ < c the relativistic velocity curves, given by Eq. I|24|l . exhibit 
essential deviations from the purely exponential decay, predicted by the Newtonian theory. 



2. Stochastic force 



We now construct a relativistic generalization of the stochastic force. To this end, we consider Eq. as an 

operational definition for the proper time parameter r. The generalization procedure will be based on the standard 
assumption (postulate) that, in temporarily comoving incrtial frames £*, the relativistic equations of motions must 
reduce to the Newtonian equations of motions. According to this assumption, for frames E», comoving with the 
particle at laboratory time t, the relativistic stochastic differential equation must reduce to 



dp*(i) = — v [p*{t) — mU*] di* + w*(t), 



(25) 



where the momentum increments w» = dlU* represent a Wiener-process with parameter D, i.e., the increments w*(t) 
have a Gaussian distribution: 



i 

V47T.D dt. 



: exp 



w*(t) 2 
TDdZ 



(26) 



1 2 3 4 5 

t [V" 1 ] 



FIG. 1: Velocity curves v(t), corresponding to Eq. 1241 . for the purely damped motion of a relativistic particle in the rest 
frame of the heat bath (laboratory frame). Especially at high velocities \v\ < c, the relativistic velocity curves deviate from the 
exponential decay, predicted by the Newtonian theory. 



Note that also in the relativistic theory the momentum increments w(t) = dW(t) may tend to infinity, as long as the 
related velocity increments remain bounded. In other words, in the relativistic theory one must carefully distinguish 
between stochastic momentum and velocity increments (this is not necessary in the nonrclativistic theory, since 
Newtonian momenta are simply proportional to their velocities). 
The next step is now to define the increment 1+1-vector by 

«) = (0,O. (27) 

This definition is in agreement with the requirement that in a comoving inertial system the 0-component of the 
1+1-force vector must vanish (see, e.g., |25| Chap. 2.3, and also compare Eqs. (|12|) . I|31l) and (|32|l of the present paper). 
Moreover, if the Lorentz frame E* is comoving with the Brownian particle at given time t, then the (equal-time) white 
noise relations I0J generalize to 

«(*)) = 0, «(*) «*(*)) = a = ° . and/ ° r " = °' (28) 

1 ' ' 1 • ! ' , otherwise. 




The rhs. of the second equation in 128fl makes it plausible to introduce a correlation tensor by 

)-f° 
~* l3 > \0 2Ddt 

thus, 

(w^(t)wUt))=D^. (29b) 
Additionally defining an 'inverse ' correlation tensor by 

(2 (2D^)"0' (29C) 
allows us to generalize the distribution of the increments from Eq. 112611 as follows 

1 



Vl +1 [w*(t)} = - 1 exp 



2 D* a/s w°(t)wUt) 



x6[w°(t)}- (30) 



Here, the Dirac (5-function on the rhs. accounts for the fact that the 0-component of the stochastic force must vanish 
in every inertial frame, comoving with the Brownian particle at time t; compare Eq. (|27|) . This also follows more 
generally from the identity 

Q=^-(-mc 2 )=m^-{u a u a ) = 2u a f a , (31) 

which, in the case of the stochastic force, translates to 

= u a w a . (32) 
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Hence, we can rewrite the probability distribution (|30(l as 

1 



Kit)] 



exp 



x 8 «;?(*)] , 



(33) 



where (it* a ) = (— c,0) is the covariant 1+1- velocity of the particle a the comoving rest frame. It should be stressed 
that, because of the constraint l|32|l . only one of the two increments w a = dW a is to be regarded as 'independent', 
which is reflected by the appearance of the 6- function in l|33|) . Also note that, due to the prefactor c, the normalization 
condition takes the simple form 



n J —CO 



Furthermore, analogous to l|17f) . we have the following more general representation of the correlation tensors 



D a/ 3 = 2D dr (ri a [_ 

6afi = ^d7 



UgU(3 

c 2 
u a up 



(34) 

(35a) 
(35b) 



Then, in an arbitrary Lorentz frame, the density (|33|) can be written as 



V 1+1 [w a (T)} = 



V 4ttD dr 

c 

dr 



exp 
exp 



- b a pW a {T)w l3 {T) 



W a {T) iii°(t) 
4L>dr 



x 6 [u a w a (T)} 
x S [u a /(t)] . 



(35c) 



To obtain the last line from the first, we have inserted D a p from (|35b(l and then used that u a w a = 0, see Eq. (|32() . 

By virtue of the above results, we are now in the position to write down the covariant Langevin equations with 
respect to an arbitrary inertial system: If a Brownian particle with rest mass m, proper time r and 1+1-velocity 
is surrounded by an isotropic, homogeneous heat bath with constant 1+1-velocity U@ , then the relativistic Langevin 
equations of motions read 



dx a (r) = ^^dr 



m 



dp a {T) = -v a p [/(t) - mU 13 ] dr + w a (r), 
where, according to Eq. i|17|) . the friction tensor is given by 

U a Ufj 



v" = v\ r p + 



(36a) 
(36b) 

(36c) 



with v denoting the viscous friction coefficient measured in the rest frame of the particle. This is a first main result of 
this work. The stochastic increments w a (r) = dF+ Q (r) are distributed according to (|35c|l and, therefore, characterized 
by 



K(r)) = 0, 



(w q (t)u/V)) = 



0, 



D a0 , T = 



(36d) 
(36e) 



with D q p given by l|35a(l . Note that in each comoving Lorentz frame, in which, at a given moment t, the particle is 
at rest, the marginal distribution of the spatial momentum increments, defined by 



d[w°(t)]V 1+l [w a (t)], 



(37) 



reduces to a Gaussian. In the Newtonian limit case, corresponding to v 2 <C c 2 , one thus recovers from Eqs. I|35|) and 
(|36ll the usual nonrelativistic Brownian motion. 
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C. Langevin dynamics in the laboratory frame 



A laboratory frame So is, by definition, an inertial system, in which the heat bath is at rest, i.e., in So we have 
(U 13 ) = (c, 0) for all times t. Hence, with respect to S -coordinates, the two stochastic differential Eqs. I|36b|l take 
the form 



dp = — vpdt + w(t), 
dE = -vpvdt + cw°(t). 



(38a) 
(38b) 



Here it is important to notice that the stochastic increments w a (t), appearing on the rhs. of l|38(l . are not of simple 
Gaussian type anymore. Instead, their distribution now also depends on the particle velocity v. This becomes 
immediately evident, when we rewrite the increment density (I35cl) in terms of So-coordinates. Using 



(u a ) = (-7C,7«) , 



we find 



V 1+1 [w a {t)]=c 



7 



AttD dt 



1/2 



exp 



w(tf - w°(t) 2 



AD df/7 



(10°) = (w°,w) 



x 5 [cjw°(t) — jvw(t)] 



(39) 



(40) 



As we already pointed out earlier, the (5-function in l|40|l reflects the fact that the energy increment w° is coupled to 
the spatial (momentum) increment w via 



= u a w a = —cyw + ryvw => w° 

Hence, w° can be eliminated from the Langevin equations l|38bfl . yielding 

dE = —vpvdt + v w(t) = v dp. 

Using the identity 



vw 
c 



cp 



\J m 2 c 2 + p 2 



we can further rewrite (|42t as 

dE = 



cp 



\J m 2 c 2 + p- 



:dp 



E{t) = v / m 2 c 4 +p{ty C 



2„2 



(41) 
(42) 
(43) 

(44) 



Thus, in the laboratory frame So the relativistic Brownian motion is completely described by the Langevin equation 
(|38a|l already. If we assume that the Brownian particle has fixed initial momentum p(0) = po or initial velocity 
v(0) = i>o, respectively, then the formal solution of l|38a(l reads (|2(J Chap. IX. 1) 



p{t) = p e 



e vs w{s). 



(45) 



The stochastic process l)45|l is determined by the marginal distribution P 1 [?/;(£)], defined in Eq. (|37|l . Performing the 
integration over the ^-function in l|4U|) . we find 



1 



47tD7 dt 



1/2 



exp 



w{tf 



where 



1 - 



-1/2 



1 



¥ 



ADi dt 

2 nl/2 



(46a) 



(46b) 



On the basis of Eqs. 138a|l and 1)46(1 one can immediately perform computer simulations, provided one still specifies 
the rules of stochastic calculus, i.e., which value of p is to be taken to determine 7 in (|46|l . In Sec. HVI several numerical 
results are presented. Before, it is useful to consider in more detail the Fokker-Planck equations of the relativistic 
Brownian motion in the laboratory frame So. By doing so in the next section, it will become clear that, for example, 
choosing p = pit) in Eqs. (|46|l . would be consistent with an Ito-interpretation |2H I33L l34l | of the stochastic differential 
equation l|38a|l . However, we will also see that alternative interpretations lead to reasonable results as well. 
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III. DERIVATION OF CORRESPONDING FOKKER-PLANCK EQUATIONS 

The objective in this part is to derive relativistic Fokker-Planck equations (FPE) for the momentum density f(t,p) 
of a free particle in the laboratory frame S . Before we deal with this problem in Sec. IIII Bl it is useful to briefly 
recall the nonrelativistic case. 



A. Nonrelativistic case 



Consider the nonrelativistic Langevin equation Ijlbjl 



(47a) 



where pit) = mv(t) denotes the nonrelativistic momentum, and, in agreement with (|3c|l . the Langevin force L(t) is 
distributed according to 



1/2 



exp 



-— Hi? 



(47b) 



As is well known |20L l35| . the related momentum probability density f(t,p) is governed by the Fokker-Planck equation 



d 



d 



dt dp 

whose stationary solution is the Maxwell distribution 







-J= — \ V V! + D — f 



m = ( 



2nD 



1/2 



up 2 \ 
eXP (-25 I 



(48) 



(49) 



B. Relativistic case 



We next discuss three different relativistic Fokker-Planck equations for the momentum density f(t,p), related to 
the stochastic processes defined by (|38a[l and (|46|) . 

Our starting point is the relativistic Langevin equation i|38a[) . which holds in the laboratory frame So (Le. in the 
rest frame of the heat bath). Next we define a new stochastic process by 



and using (|46 bj l . we can rewrite l|38a[l as 

dp = -vpdt + y/j y(t), 
where y(t) is distributed according to the momentum- independent density 



( i 



\4ttD dt / 



v 1/2 


[ v(tf] 


J exp 


ADdt 



(50) 



(51a) 



(51b) 



Thus, instead of the increments w(t), which implicitly depend on the stochastic process p via Eqs. I|4ti|) . we consider 
ordinary p- independent white noise y(t), determined by l|51b|) . from now on. Due to the multiplicative coupling of 
y(t) in (|51a|l . we must next specify rules for the "multiplication with white noise" [note that, upon viewing Eqs. Ullfl . 
(13511 and (|36|l as postulates of the relativistic Brownian motion, all above considerations remain valid independent of 
this specification]. 

In the followin g su bsections, we shall discuss thr ee popu lar multiplication rules, which go back to proposals made 
by Ito ElllS mill, by Stratonovich and Fisk [H HI El H HI , and by Hanggi HHHgl and Klim ontovich [p, 
respectively. As well-known [2(], [3(j , these different interpretations of the stochastic process (|51f) result in different 
Fokker-Planck equations, i.e., the Langevin equation (|51|l per se does not uniquely determine the corresponding 
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Fokker-Planck equation; it is the stochastic interpretation of the multiplicative noise that matters from a physical 
point of view. 

Nevertheless, the three approaches discussed below have in common that, formally, the related Fokker-Planck 
equation can be written as a continuity equation |4^ : 

p(t,p) + ^j(t,p)=0, (52) 

but with different expressions for the probability current j(t,p). It is worthwhile to anticipate that only for the 
Hanggi-Klimontovich approach (see Sec. IIII B 3|) the current j(t,p) takes such a form that the stationary distribution 
of l|52l) can be identified with Jiittner's relativistic Maxwell distribution [24j. 



1. Ito approach 



According to Ito's interpretation of the Langevin equation l|51a|) . the coefficient before y(t) is to be evaluated at 
the lower boundary of the interval [t, t + dt], i.e., we use the pre-point discretization rule 



7 = l{p{t)); 



(53) 



where as before 



lip) = 1 



1/2 



Ito's choice leads to the following expression for the current [19l.l20l.l33l 13 



3i(p,t) = - 



d 

vpf + D—-y(p)f 



(54) 



The related relativistic Fokker-Planck equation is obtained by inserting this current into the conservation law H52fl . 
The current 154fl vanishes identically for 



lip) 



-5/* 



lip) 



(55) 



where C\ is the normalization constant. Consequently, fi{p) is a stationary solution of the Fokker-Planck equation. 
In view of the fact that 



dp . — c 2 m 2 \ 1 1 
liP) 



we find the following explicit representation of i|55|) 



V 



(56) 



hip) =C I (l + 



-1/2 



exp -0\ 1 + 



(57) 



where 







9 9 
vm c 

D 



(58) 



The dimcnsionless parameter (3 can be used to define the scalar temperature T of the heat bath via the Einstein 
relation 



nc D 

/3 vnv ' 



(59) 



with ks denoting the Boltzmann constant. Put differently, the parameter = mc 2 /(fcgT) measures the ratio between 
rest mass and thermal energy of the Brownian particle. 
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2. Stratonovich approach 



According to Stratonovich, the coefficient before y{t) in (|51a|l is to be evaluated with the mid-point discretization 
rule, i.e., 



7 = 7 



/ P (t) + P (t + at) 



This choice leads to a different expression for the current [1 

js(p,t) = - 



namely 



i<P/ + £>Vt^)— VtW / 



This Stratonovich-Fisk current js vanishes identically for 

C s 



/s(p) 



vt(p) 



exp 



l)J dp W) 



and, by virtue of l|56|l . the explicit stationary solution of Stratonovich's Fokker-Planck equation reads 

-1/4 / 



fs(p) = C s 1 + 



P 



exp -/3W1 + 



P 



(60) 



(61) 



(62) 



(63) 



3. Hdnggi-Klimontovich approach 

Now let us still consider the Hanggi-Klimontovich stochastic integral interpretation, sometimes referred to as the 
transport form [ifl l4ll or also as the kinetic form [3lJ ■ According to this interpretation, the coefficient in front 
of y{t) in (|51a|l is to be evaluated at the upper boundary value of the interval [t, t+ di\; i.e., within the post-point 
discretization we set 



7 = 7 (p(i + di)). 

This choice leads to the following expression for the current [3ll l4ll E 



d ' 

vpf + Dj(p)—f 



Jhk(p,*) = 
The current j'hk vanishes identically for 

/hk(p) = Cmcexp 
and, by virtue of H56[). the stationary solution explicitly reads 

/hk(p) = C*hk exp -j3\/l 



17 / ^ p 



(64) 



(65) 



(66) 



(67a) 



Using the temperature definition in l|59|) and the relativistic kinetic energy formula E = \J m 2 ^ + p 2 c 2 , one can 
further rewrite l|67a() in a more concise form as 



/hk(p) = Chkgxp - 



k B T 



(67b) 



The distribution function l|67l) is known as the relativistic Maxwell distribution. It was first obtained by F. Jiittner [24j 
back in 1911. Pursuing a completely different line of reasoning, he found that l|67|l describes the velocity distribution of 
the non- interacting relativistic gas (see also [43^- In contrast to our approach, which started out with constructing the 
relativistic generalization of the Langevin equations, Juttner's derivation started from a maximum-entropy-principle 
for the gas. 
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By comparing i|55|) . I|63|) and i|67a[) one readily observes that the stationary solutions fi/s differ from the Jiittner 
function /hk through additional p-dependent prefactors. In order to illustrate the differences between the different 
stationary solutions, it useful to consider the related velocity probability density functions (/>i/s/hk(w), which can be 
obtained by applying the general transformation law 



(v) = f(p(v)) 



dp 
dv 



(68) 



in combination with 



mv 



P = 



y/1 - V 2 / C 2 ' 



The exponential factor ensures that the velocity density functions ^i/s/hk^) are in fact zero if v 2 > c 2 . 

In Fig. |2]we have plotted the probability density functions ^i/s/hk(") f° r different values of the parameter (3. The 
normalization constants were determined by numerically integrating 4>{v) over the interval [— c, c]. As one can observe 
in diagram [3 (a), for large values of /3, corresponding to small temperature values fcgT <C mc 2 , the density functions 
•^i/S/hk^) approach a common Gaussian shape. On the other hand, for high temperature values kgT > mc 2 the 
deviations from the Gaussian shape become essential. The reason is that, for a (virtual) Brownian ensemble in the 
high-temperature regime, the majority of particles assumes velocities, that are close to the speed of light. It is also 
clear that in other Lorentz frames X', which are not rest frames of the heat bath, the stationary distributions will no 
longer stay symmetric around v = 0. Instead, they will be centered around the non-vanishing E'-velocity V of the 
heat bath. 
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FIG. 2: Stationary solutions 4>j/s/-rk{ v ) of the relativistic Fokker-Planck equations, according to Ito (hsolid line), Stratonovich 
(S:dotted) and Hanggi-Klimontovich (HKdashed-dotted). For low temperatures, i.e., for |J > 1, a Gaussian shape is ap- 
proached, see diagram (a). On the other hand, for very high temperature values, corresponding to f3 < 1, the distributions 
exhibit a bistable shape, and the quantitative deviations between <^i/s/hk(") increase significantly as f3 — » 0. 



As an obvious question then arises, which of the above approaches (Ito, Stratonovich or Hanggi-Klimontovich) is 
the physically correct one. We believe that, at this level of analysis, it is impossible to provide a definite answer to this 
question. Most likely, the answer to this problem requires additional information about the microscopic structure of 
the heat bath (see, e.g., the discussion of Ito-Stratonovich dilemma in the context of 'internal/external' noise as given 
in Chap. IX. 5 of van Kampen's textbook [2(|). At this point, it might be worthwhile to mention that the relativistic 
Maxwell distribution (|BT|) is also obtained via the transfer probability method used by Schay, see Eq. (3.63) and 
(3.64) in Ref. 0, and that this distribution also results in the relativistic kinetic theory [2|j. By contrast, the recent 
work of Franchi and Le Jan is based on the Stratonovich approach. From pure physical insight, however, it is the 
transport form interpretation of Hanggi and Klimontovich that is expected to provide the physical correct description. 



IV. NUMERICAL INVESTIGATIONS 



The numerical results presented in this section were obtained on the basis of the relativistic Langevin equation 
(|51|l . which holds in the laboratory frame £o- F° r simplicity, we confined ourselves here to considering the Ito- 
discretization scheme with fixed time step dt , see Sec. IIII B II In all simulation we have used an ensemble size of 
N = 10000 particles. Moreover, a characteristic unit system was fixed by setting m = c = v = 1. Formally, this 
corresponds to using re-scaled dimensionless quantities, such as p = p/mc, x = xv/c, t — tv, v — v/c etc.. The 
simulation time-step was always chosen as dt = 0.001^ -1 , and the Gaussian random variables y(t) were generated by 
using a standard random number generator. 
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A. Distribution functions 



In our simulations we have numerically measured the cumulative velocity distribution function F(t, v) in the lab- 
oratory frame Eg. Given the probability density <p(t,v), the cumulative velocity distribution function is defined 
by 

/V 
du(f>(t,u). (69) 

In order to obtain F(t, v) from numerical simulations, one simply measures the relative fraction of particles with 
velocities in the interval [— c, v). Figure[31shows the numerically determined stationary distribution functions (squares), 
taken at time t = lOOi^ 1 and also the corresponding analytical curves -Fi/s/hk^)- The latter were obtained by 
numerically integrating the formula (|69J) . using the three different stationary density functions </>i/s/hk(^) from Sec. 

eh 




FIG. 3: These diagrams show a comparison between numerical and analytical results for the stationary cumulative distribution 
function F(v) in the laboratory frame £q. (a) In the nonrelativistic limit /3 > 1 the stationary solutions of the three different 
FPE are nearly indistinguishable, (b-c) In the relativistic limit case < 1, however, the stationary solutions exhibit deviations 
from each other. Because our simulations are based on an Ito-discretization scheme, the numerical values (squares) are best 
fitted by the Ito solution (solid line). 



As one can see in diagram|21(a), for low temperature values, corresponding to (3 ^> 1, the three stationary distribution 
functions are nearly indistinguishable. For high temperatures, corresponding to [3 < 1, the stationary solutions exhibit 
significant quantitative differences, see Fig. 03(b) and (c). Since our simulations are based on an Ito-discretization 
scheme, the numerical values (squares) are best fitted by the Ito solution (solid line). Also note that the quality of 
the fit is very good for the parameters chosen in the simulations, and that this property is conserved over several 
magnitudes of (3. This suggests that numerical simulations of the Langevin equations provide a very useful tool, if 
one wishes to study relativistic Brownian motions in more complicated settings, e.g., in higher dimensions or in the 
presence of additional external fields and interactions. In this context, it should again be stressed that the appropriate 
choice of the discretization rule is especially important in applications to realistic systems. 



B. Mean square displacement 



In this part we consider the spatial mean square displacement of the free relativistic Brownian motion. Since this 
quantity is easily accessible in experiments, it has played an important role in the verification of the nonrelativistic 
theory. 

As before, we consider an ensemble of N independent Brownian particles with coordinates x^(t) in So and initial 
conditions Xu\ (0) = 0, v^) (0) = for i = 1, 2, . . . , N. The position mean value is defined as 

1 N 

^ = ^1^)^ ( 7 °) 

i=l 

and the related second moment is given by 

— 1 N 
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The empirical mean square displacement can then be defined as follows 

a 2 (t)^(t)-[x(t)} 2 . 

Cornerstone results in the nonrelativistic theory of the one-dimensional Brownian motion are 

lim x(t) — ► 0, 

t — >+oo 

a 2 (t) 



lim 

t — > ~\~ OO 



t 



2D X , 



where the constant 



D x = 



D 



(72) 

(73a) 
(73b) 

(74) 



is the nonrelativistic coefficient of diffusion in coordinate space (not to be confused with noise parameter D). 

It is therefore interesting to consider the asymptotic behavior of the quantity a 2 (t) / 1 for relativistic Brownian mo- 
tions, using again the Ito-relativistic Langevin dynamics from Sec. IIII B II In Fig. 01(a) one can see the corresponding 
numerical results for different values of j3. As one can observe in this diagram, for each value of (3, the quantity a 2 (t) /t 
converges to a constant value. This means that, at least in the laboratory frame So, the asymptotic mean square 
displacement of the free relativistic Brownian motions increases linearly with t. For completeness, we mention that 
according to our simulations the asymptotic relation l|73a|l holds in the relativistic case, too. 

In spite of these similarities between nonrelativistic and relativistic theory, an essential difference consists in the 
explicit temperature dependence of the limit value 2D X . As illustrated in Fig. El(b), the numerical limit values 2Df Q0 , 
measured at time t = lOOv^ 1 , arc well fitted by the empirical formula 



(75) 



1/09 + 2)' 

which reduces to the nonrelativistic result (|74|l in the limit case j3 ^> 2 (low-temperature limit case). We will leave it 
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FIG. 4: (a) Mean square displacement, divided by time t, as numerically calculated for different /3- values in the laboratory frame 
£0 (rest frame of the heat bath). As evident from this diagram, for the relativistic Brownian motion the related asymptotic mean 
square displacement grows linearly with t. (b) The coordinate space diffusion constant D x 00 (f3) was numerically determined at 
time t = \00v~ 1 . The dashed line corresponds to the empirical fitting formula D x ((3) = <?v~ x (fi + 2) _1 , which reduces to the 



classical nonrelativistic result D x 



y(vfi) = kT/{mv) for /3 > 2. 



as an open problem here, to find an analytical justification for the empirically determined formula (|75[) . Instead we 
merely mention that, upon noting 143(1 . the relativistic Fokker-Planck equations for the full phase space density reads 

m f(f ' P ' x) + ^ m Z + p 2 i f{t,P,X) = - j^/s/HK(*,P,z), (76) 

which might serve as a suitable starting point for such an analysis. Compared with the relativistic Fokker-Planck 
equations from the preceding section, the second term on the lhs. of l(76|l is new. In particular, we recover the 
relativistic Fokker-Planck equations for the marginal density f(t,p), see Sec. IIIII by integrating Eq. (|76|l over a 
spatial volume with appropriate boundary conditions. Finally, we mention once again that also (|76|l . as well as all 
the other results that have been presented in this part, exclusively refer to the laboratory frame So- 
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V. CONCLUSION 

Concentrating on the simplest case of 1+1-dimensions, we have put forward the Langevin dynamics for the stochastic 
motion of free relativistic Brownian particles in a viscous medium (heat bath). Analogous to the nonrelativistic 
Ornstein-Uhlenbeck theory of Brownian motion [rH IHI I20I |44| . it was assumed that the heat bath can, in good 
approximation, be regarded as homogenous. Based on this assumption, a covariant generalization of the Langevin 
equations has been constructed in Sec. [n] According to these generalized stochastic differential equations, the viscous 
friction between Brownian particle and heat bath is modeled by a friction tensor v a a. For a homogeneous heat bath 
this friction tensor has the same structure as the pressure tensor of a perfect fluid [25|]. In particular, it is uniquely 
determined by the value of the (scalar) viscous friction coefficient v, measured in the instantaneous rest frame of 
the particle (Sec. Ill B 1(1 . Similarly, the amplitude of the stochastic force is also governed by a single parameter D, 
specifying the Gaussian fluctuations of the heat bath, as seen in the instantaneous rest frame of the particle (Sec. 
IT1B21) . 

In Sec. Ill CI the relativistic Langevin equations have been derived in special laboratory coordinates, corresponding 
to a specific class of Lorentz frames, in which the heat bath is assumed to be at rest (at all times). One finds 
that the corresponding relativistic distribution of the momentum increments now also depends on the history of the 
momentum coordinate. This fact is in contrast with the properties of ordinary Wiener processes 0, |2^| , underlying 
nonrelativistic standard Brownian motions with 'additive' Gaussian white noise. However, as shown in Sec. II I II it is 
possible to find an equivalent Langevin equation, containing 'multiplicative' Gaussian white noise. 

In order to achieve a more complete picture of the relativistic Brownian motion, the corresponding relativistic 
Fokker-Planck equations (FPE) have been discussed in Sec. IIIII (again with respect to the laboratory coordinates with 
the heat bath at rest). Analogous to nonrelativistic processes with 'multiplicative' noise, one can opt for different 
interpretations of the stochastic differential equation, which result in different FPE. In this paper, we concentrated on 
the three most popular cases, namely the Ito, the Stratonovich-Fisk, and the Hanggi-Klimontovich interpretation. We 
discussed and compared the corresponding stationary solutions for a free Brownian particle. It could be established 
that only the Hanggi-Klimontovich interpretation is consistent with the relativistic Maxwell distribution. This very 
distribution was first derived by Juttner [24j as the equilibrium velocity distribution of the relativistic ideal gas. Later 
on, it was also discussed by Schayin the context of relativistic diffusions Q and by de Groot et al. in the framework 
of the relativistic kinetic theory [23 . 

In Sec. IIVI we presented numerical results, obtained on the basis of an Ito pre-point discretization rule. The 
simulations indicate that - analogous to the nonrelativistic case - the relativistic mean square displacement grows 
linearly with the laboratory coordinate time; the temperature dependence of the related spatial diffusion constant, 
however, becomes more intricate. In principle, the numerical results suggest that simulations of the Langevin equations 
may provide a very useful tool for studying the dynamics of relativistic Brownian particles. In this context it has 
to be stressed that an appropriate choice of the discretization rule is esp ecially important in applications to realistic 
physical systems. If, for example, agreement with the kinetic theory [29j is desirable, then a post-point discretization 
rule should be used. 

From the methodical point of view, the systematic relativistic Langevin approach of the present paper differs from 
Schay's transition probability approach Q and also from the techniques applied by other authors 0, IE • As we 
shall discuss in a forthcoming contribution, the above approach can easily be generalized to settings which are more 
relevant with regard to experiments (such as the 1+3-dimensional case, the presence of additional external force fields, 
etc.). 

With regard to future work, several challenges remain to be solved. For example, one should try to derive an 
analytic expression for the temperature dependence of the spatial diffusion constant. A suitable starting point for 
such studies might be the FPE for the full phase space density given in Eq. (|76|l . Another possible task consists in 
finding explicit exact or at least approximate time-dependent solutions of the relativistic FPE. Furthermore, it seems 
also interesting to consider extensions to general relativity, as, to some extent, recently discussed in the mathematical 
literature |l3| . In this context, the physical consequences of the different interpretations (Ito vs. Stratonovich vs. 
Hanggi-Klimontovich) become particularly interesting. 
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